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THE  RESPONSE  OF  THE  ZONALLY 

AVERAGED  CIRCULATION  TO  STRATOSPHERIC  OZONE  REDUCTIONS 

I.  INTRODUCTION 

Considerable  effort  has  been  expended  to  access  the  impact  of  chemical  pollutants  on  the 
stratospheric  environment.  These  pollutants  catalytically  attack  ozone,  the  principal  atmos- 
pheric absorber  of  biologically  harmful  solar  UV  radiation.  This  UV  radiation  also  heats  the 
stratosphere  and  drives  the  zonally  averaged  (mean)  circulation  of  the  middle  atmosphere. 
Since  the  stratosphere  is  a relatively  stagnant  region  of  the  atmosphere,  environmental  assess- 
ment of  pollution  effects  will  depend  critically  on  the  response  of  stratospheric  dynamics  to  the 
changes  In  radiative  forcing  that  accompany  pollution  induced  ozone  reductions. 

Our  recently  developed  numerical  model  (Schoeberl  and  Strobel,  1977;  hereafter  referred 
to  as  Paper  I)  has  been  used  to  compute  the  response  of  the  mean  circulation  to  ozone  density 
perturbations  in  the  middle  atmosphere.  We  have  designed  a series  of  numerical  experiments 
to  simulate  observed  natural  perturbations  and  predicted  man-made  perturbations.  Both  the 
immediate  time  dependent  response  and  the  final  asymptotic  state  (steady  state  solution)  for 
these  perturbations  have  been  examined. 

II.  MODEL 

In  Paper  I we  previously  studied  the  zonally  averaged  circulation  of  the  middle  atmos- 
phere (15-125  km)  with  a quasigeoslrophic,  numerical  model  that  explicitly  includes  a self- 
consistent  calculation  of  solar  radiative  heating  due  to  O2  and  O3  absorption,  Newtonian  cool- 
ing, Rayleigh  friction,  and  tropopause  boundary  conditions  based  on  climatological  averages. 
Manuscript  submitted:  September  29,  1977. 
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In  order  to  calculate  the  effects  of  ozone  reduction  we  must  also  compute  the  globally  averaged 
atmospheric  fields,  denoted  by  < >,  as  well  as  (he  zonally  averaged  fields,  denoted  by  an  over- 
bar, which  are  equal  to  the  difference  between  zonally  averaged,  variables  and  globally  averaged 
variables.  ' 


The  globally  averaged  atmosphere  is  assumed  to  be  motionless,  hydrostatic,  and  in  ap- 


, where  p„  is  a reference 


proximate  radiative  balance.  With  the  vertical  coordinate  r = //;  | 
pressure  taken  to  be  100  mb,  the  governing  equations  of  the  globally  averaged  atmosphere  are 
b<T> 
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where  a is  the  Newtonian  cooling  coefficient;  b is  the  mean  cooling  rate,  H is  the  heating  rate 
per  unit  mass,  7"  is  temperature,  is  geopotential  height,  and  R is  the  dry  air  gas  constant. 

The  zonally  averaged  circulation  is  governed  by 
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where  is  the  Rayleigh  friction  coefficient,  y = Ay/211,  11  is  the  earth's  rotational  frequency, 
a is  the  earth’s  radius,  S is  the  static  stability,  and  K is  latitude.  The  mean  zonal  wind  is  given 
by 


—sin  b 


b<f> 


2ila(y^  -f  sin^w)  6w 

The  Newtonian  cooling  coefficient,  <<«>,  is  based  on  Dickinson’s  (1973)  detailed  IR  cal- 
culations and  is  assumed  invariant  on  pressure  surfaces.  Since  (he  globally  averaged  tempera- 
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ture  profile  is  fairly  well  known  (CIRA,  1972)  and  the  major  contributions  to  solar  heating, 
<//>,  can  be  accurately  and  self-consistently  computed  as  described  in  Paper  I,  <b>  is 

d ^ Z'  ^ 

chosen  to  obtain  radiative  balance  for  the  CIRA  (1972)  <T>  profile  when  r = 0 in 

at 

Eq.  (1).  Below  70  km  our  derived  <b>  values  are  in  accord  with  the  Dickinson  (1973) 
results.  These  < b>  values  are  also  assumed  invariant  (with  time  and  changes  in  atmospheric 
structure)  on  pressure  surfaces  in  all  calculations.  The  boundary  conditions  used  for  the  nu- 
merical solution  of  Eq.  (3)  are  given  in  Paper  I.  For  small  O3  density  perturbations,  it  is  a rea- 
sonable assumption  that  the  IR  radiative  transfer  and  hence  cooling  rates  will  not  be  seriously 
altered.  But  for  large  O3  density  perturbations  our  computations  should  be  regarded  as  only 
qualitatively  correct.  The  photochemical  acceleration  of  the  thermal  relaxation  rate  was 
neglected  in  Eqs.  (1)  and  (3)  because  the  densities  were  specified,  the  catalytic  reactions 
have  a weak  temperature  dependence,  and  the  rate  depends  on  the  nature  of  the  chemical  per- 
turbation (pollutant).  The  neglect  of  this  rate  will  only  affect  the  time  dependent  response  of 
the  stratopause  region. 

III.  UNIFORM  OZONE  DENSITY  REDUCTIONS 

Probably  the  simplest  experiment  which  can  be  performed  is  to  uniformly  reduce  ozone 
everywhere  Three  such  experiments  were  performed  with  a 25%,  50%,  and  75%  ozone  reduc- 
tion, respectively.  Selected  asymptotic  values  of  the  meteorological  variables  are  presented  in 
Table  1 as  diagnostic  indicators  of  the  mean  circulation. 

With  any  uniform  ozone  reduction  two  important  effects  occur:  first,  the  globally  aver- 
aged stratospheric  temperature  decreases  which  causes  pressure  surfaces  to  descend  in  altitude. 
Second,  the  differential  heating  rate  that  drives  the  mean  circulation  decreases  which  leads  to  a 
smaller  meridional  temperature  gradient  and  a weaker  polar  night  jet. 
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These  effects  are  quantitatively  displayed  in  Figs.  1 and  2 where  the  globally  averaged 


heating  and  temperature  profiles,  respectively,  are  preser  The  pole  to  pole  diffe  ;ntial  heat- 
ing rate  is  proportional  to  the  globally  averaged  heating  rate,  <H>.  We  note  in  Fig.  2 that  the 
stratopause  and  mesopause  descend  significantly  as  O3  density  reductions  increase.  Also,  radia- 
tion in  the  Hartley  and  Huggins  bands  penetrates  deeper  and  heats  a denser  region  of  the  at- 
mosphere. At  unit  optical  depth  the  peak  heating  rate  per  unit  mass  is  smaller,  as  illustrated  in 
Fig.  1.  The  globally  averaged  temperature  profiles  clearly  reflect  this  effect  in  Fig.  2.  A lower 
mesosphere  results  from  the  smaller  vertical  extent  of  a cooler  stratosphere. 

Figure  3 shows  the  structure  of  the  polar  night  jet  for  each  of  the  three  experiments 
along  with  the  unperturbed  atmosphere.  Note  that  the  height  coordinates,  which  represent  the 
globally  averaged  height  of  individual  pressure  surfaces,  are  different  for  each  figure  as  a conse- 
quence of  the  altitude  variation  of  pressure  surfaces. 

Ozone  reduction  leads  to  not  only  a cooler,  collapsing  stratosphere,  but  by  the  thermal 
wind  relation,  to  descending  zonal  jets  which  decrease  in  intensity  as  seen  in  Figure  3.  Note 
that  as  the  jets  weaken  the  influence  of  tropospheric  height  fields  extends  slightly  higher  into 
the  stratosphere.  This  is  especially  apparent  in  the  summer  hemisphere. 

From  steady  state  calculations  the  stratospheric  response  time  to  perturbations  can  not  be 
inferred,  although  the  damping  rates,  and  <«>,  would  be  expected  to  yield  a characteristic 
relaxation  time.  A time  dependent  calculation  with  a 75"/!  O3  density  reduction  switched  on  at 
f “ 0 is  illustrated  by  the  changes  in  i7  at  62.5°N  and  Tat  87.5‘’N  (winter  hemisphere)  at  vari- 
ous constant  altitudes  in  Fig.  4.  The  e-folding  rate  for  the  zonal  wind  (Fig.  4a'  is  ~5  days  in 
the  mesosphere  and  upper  stratosphere  and  somewhat  longer  in  the  middle  stratosphere.  The 
temperature  perturbations  in  Fig.  4b  are  more  complicated  as  a result  of  the  collapsing  stratos- 
phere. For  example  at  95  km  the  temperature  initially  decreases  from  the  reduced  heating,  but 
eventually  increases  as  the  warm  thermosphere  descends  to  lower  levels.  Similarly  the  warm 
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winter  mesosphere  descends  and  negates  a large  temperature  decrease  in  the  vicinity  of  the 


stratopause  at  55  km.  At  this  altitude  a relaxation  time  of  —5  days  could  be  deduced,  although 
admittedly  the  interference  of  lime  constants  from  Fig,  4b  is  not  straightforward.  Another  in- 
teresting feature  of  Fig.  4 is  the  slight  increase  in  i7and  Fafler  15  days.  This  oscillation  in  the 
response  is  a result  of  the  interaction  of  the  globally  averaged  temperature  with  the  zonal  mean 
perturbations  through  the  static  stability,  S.  As  < T>  decreases,  5' also  decreases  and  the  zonal- 
ly  averaged  temperature  approaches  radiative  equilibrium.  As  shown  in  Paper  1,  the  polar  night 
jets  are  much  stronger  for  a radiative  equilibrium  atmosphere.  Thus  this  effect  tends  to  coun- 
teract the  decrease  in  the  meridional  heating  gradient  which  weakens  the  pole  to  pole  circula- 
tion. 

IV.  REALISTIC  OZONE  REDUCTIONS 

The  results  of  the  previous  computer  experiments  can  not  be  realistically  applied  to  the 
stratosphere  since  the  column  O3  density  reductions  are  substantially  larger  than  both  ob- 
served and  predicted  perturbations  of  the  ozone  layer  (Heath  et  al.,  1977;  Committee  on  Im- 
pact on  Stratospheric  Change  (CISC),  1976).  They  do,  however,  serve  as  benchmarks  to  evalu- 
ate pollution  effects.  To  provide  a realistic  assessment  of  the  potential  impact  on  stratospheric 
dynamics  two  additional  computer  experiments  were  performed.  Heath  et  al.  (1977)  observed 
ozone  reductions  during  the  August  1972  solar  proton  event.  A maximum  O3  density  reduc- 
tion of  —22%  was  observed  at  44  km  and  confined  to  polar  regions.  A suitable  model  based  on 
their  data  was  constructed. 

A time  dependent  simulation  of  the  dynamical  response  of  the  middle  atmosphere  to  the 
solar  proton  event  showed  negligible  impact.  The  strength  of  the  polar  night  Jet  remained  con- 
stant (less  than  1%  change)  and  the  globally  averaged  stratospheric  temperature  changed  by 
less  than  IK  at  45  km.  In  the  polar  regions  temperature  changes  of  at  most  2K  are  calculated. 
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which  is  consistent  with  the  crude  estimates  of  Zerefos  and  Crulzen  (1975).  This  result  is  not 
unexpected  in  light  of  the  results  from  the  previous  section;  the  Oj  reduction  is  confined  to 
polar  regions  and  globally  the  chemical  ozone  destruction  is  negligible. 

The  second  experiment  was  designed  to  represent  the  scenarios  predicted  by  CISC  (1976) 
in  their  evaluation  of  halocarbon  impact  on  the  ozone  layer.  The  "no  chlorine  nitrate  forma- 
tion" curve  in  their  Fig.  8.6  was  approximated  by  a suitable  linear  fit.  They  predicted  a peak 
global  ozone  reduction  of  50%  at  40  km. 

The  halocarbon  induced  O3  reduction  is  projected  to  occur  slowly  over  decades  and  a 
steady  state  treatment  should  suffice.  The  steady  state  results  are  indicated  in  Figs.  1-3  by  the 
curves  labelled  "halocarbon".  Substantial  changes  in  the  globally  averaged  heating  rate  and 
temperature  are  evident  below  55  km  where  O3  reduction  becomes  significant.  The  maximum 
perturbation  in  < T>  is  13K  at  43  km.  The  overall  jet  structure  is  similar  to  the  uniform  25% 
reduction  case,  but  the  amplitude  change  is  not  as  great. 

Even  though  the  column  O3  density  reduction  in  both  cases  is  25%,  the  column  O3  den- 
sity reduction  for  the  halocarbon  case  is  much  less  than  25%  above  48  km  where  the  peak 
heating  occurs  due  to  the  O3  Hartley  lands.  This  effect  produces  the  major  differences  between 
the  <H>  and  < 7>  profiles  for  the  25%  reduction  and  halocarbon  experiments.  Since  the 
25%  case  reaches  the  same  maximum  heating  rate  at  a lower  altitude,  the  polar  night  and  day 
jets  have  slightly  lower  intensity  because  the  vertical  shear  generated  from  thermal  wind  rela- 
tion is  sustained  over  a shorter  height  interval. 

V.  SOLAR  FLUX  CHANGES 

In  the  same  vein  as  the  uniform  ozone  density  reductions  we  have  examined  the  ..lodel 
sensitivity  to  changes  in  the  solar  flux.  Figure  5 shows  the  response  of  various  model  variables 
to  altered  solar  output.  This  calculation  ignores  the  probably  weaker  adjustment  of  O3  concen- 
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traiions;*  however,  we  believe  that  is  representative  of  the  effects  of  small  variations  in  solar 
flux. 

It  is  apparent  that  the  response  of  the  zonally  averaged  circulation  is  nearly  linear:  both 
the  globally  averaged  straiopause  temperature  and  the  strength  of  the  mean  zonal  winds  vary 
uniformly  with  flux,  At  lower  altitudes  the  response  is  weaker  as  the  tropopause  boundary  con- 
ditions are  held  fixed  for  these  calculations. 

The  amplitudes  of  planetary  waves  are  also  shown  in  Table  1 for  the  calculated  zonal  jet 
structure  variation  with  solar  flux.  For  m = 1,  the  amplitude  varies  nearly  linearly  with  jet 
strength  but  its  phase  is  insensitive  to  jet  strength  with  at  most  a 7°  shift.  This  is  roughly  the 
same  magnitude  as  the  phase  change  between  pressure  levels  of  the  model.  The  response  of 
m = 2 planetary  wave  is  not  linear  with  jet  strength  and  has  a maximum  amplitude  for  the  86 
ms  jet  (cf.  Table  1).  Changes  in  the  solar  flux  on  the  order  of  1%  would  thus  negligibly 
alter  the  large  scale  wave  structure  in  the  upper  stratosphere  and  have  an  even  smaller  effect  at 
lower  altitudes. 

VI.  CONCLUDING  REMARKS 

We  have  used  a numerical  quasigeostrophic  model  of  the  mean  circulation  in  the  middle 
atmosphere  to  investigate  the  dynamical  response  to  reductions  in  ozone  density.  The  August 
1972  solar  proton  event  produced  measurable  reductions  in  the  polar  ozone  content  but  minor 
ozone  destruction  on  a global  basis.  Our  calculations  showed  that  such  an  event  produces 
negligible  changes  in  the  mean  circulation.  Halocarbons  are  predicted  to  have  a global  impact 
on  the  ozone  layer  (CISC  , 1976)  and  we  find  approximately  10%  reduction  in  the  zonal  jet 
strength  with  smaller  percentage  changes  in  the  globally  averaged  temperature.  For  these  real- 
istic ozone  reductions  the  Dickinson  (1973)  Newtonian  cooling  coefficients  should  adequately 
represent  the  IR  radiative  transfer. 

*J  Chang  (pnvaie  tomniunicaiion) 
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To  produce  catastrophic  effects  on  the  mean  circulation  we  investigated  arbitrary  uniform 
ozone  reductions.  For  50%  and  75%  ozone  reductions  pronounced  changes  in  the  zonal  jet 
structure  and  globally  averaged  temperature  are  obtained.  It  is  unlikely  that  the  results  are 
quantitatively  accurate  since  the  IR  cooling  representation  probably  breaks  down  for  these 
large  perturbations.  We  have  also  ignored  explicit  eddy  transport  effects  by  planetary  waves  in 
the  middle  atmosphere,  since  they  are  not  yet  well  understood. 

Our  results  cannot  be  interpreted  as  a purely  linear  response  of  the  atmosphere  to  a 
reduction  in  the  differential  heating.  Changes  in  the  globally  averaged  thermal  structure  alter 

the  static  stability  which  is  nonlinearly  coupled  to  the  vertical  velocity.  In  Fig.  6 we  illustrate 
the  effect  on  the  globally  averaged  static  stability  for  the  various  ozone  reduction  cases  con- 
sidered. It  should  be  remembered  that  the  mean  circulation  has  a vertical  scale  of  approximate- 
ly three  scale  heights  (Leovy,  1964)  and  will  not  respond  to  "small  scale"  structure  in  the  static 
stability.  As  a consequence  the  non-linear  coupling  is  somewhat  less  than  one  might  infer 
from  Fig.  6. 

Hines  (1974)  suggested  that  solar-weather  interactions  may  occur  by  changes  in  the 
reflection  and  transmission  properties  of  the  middle  atmosphere  for  planetary  waves  with  solar 
activity.  The  August  1972  solar  proton  event  would  certainly  be  an  extreme  case  of  solar  ac- 
tivity affecting  the  middle  atmosphere.  Although  local  ozone  reductions  of  ~22“/i  occurred  in 
polar  regions,  the  net  effect  on  the  mean  circulation  was  negligible.  To  confirm  that  the  effect 
on  the  planetary  wave  structure  was  also  negligible  we  performed  computations  similar  to 
Schoeberl  and  Geller  (1977)  for  the  perturbed  and  unperturbed  atmospheres.  The  maximum 
stratospheric  change  in  either  amplitude  or  phase  was  less  than  0.5')(i.  In  addition  we  note  that 
the  response  time  of  the  mean  circulation  is  probably  too  long  to  explain  solar-weather  cou- 
pling, since  the  response  time  has  the  same  time  scale  as  the  normal  transient  variation  of 
planetary  waves  in  the  lower  stratosphere  (Hirota  and  Sato,  1969),  A similar  calculation  was 
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performed  for  the  halocarbon  case  with  resultant  change  in  wave  amplitude  of  about  15%  at  50 
km  for  w =1  and  10%  for  m =2  at  the  same  altitude.  This  deviation  is  still  less  than  the  nor- 
mal winter  fluctuations  of  planetary  wave  amplitude.  We  many  conclude  from  these  calcula- 
tions that  the  transport  and  mixing  of  minor  constituents  in  the  stratosphere  which  occurs  pri- 
marily through  planetary  scale  eddies  will  not  be  substantially  altered  by  predicted  ozone 
reductions. 
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Fig.  1 - The  globally  averaged  heating  rate  in  K day'^ 
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Fig.  2 - The  globally  averaged  temperature  in  K 
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LATITUDE 
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Fig.  3 - The  steady  state  mean  zonal  wide  structure  of  the  upper 
atmosphere  in  ms‘  . Part  a is  the  “normal”  structure.  Parts  h,  c, 
and  d show  the  structure  for  25%,  50%,  and  75%  reduction  in 
the  ozone  mixing  ratio,  respectively.  Part  e shows  the  results 
for  the  halocarbon  case. 
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Fig.  3 (Continued)  - Ttie  steady  state  mean  zonal  wide  structure  of  the  upper  atmosphere  in  ms"* . Part  a is  the 
“normal"  structure.  Parts  b,  c,  and  d show  the  structure  for  25%,  50%,  and  75%  reduction  in  the  ozone  mixing  ratio, 
respectively.  Part  e shows  the  results  for  the  halocarbon  case. 
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I ig.  5 - The  response  of  the  stratosphere  and  mesosphere  to  solar  flux  changes.  <T>  is  the  stratopause  temperature, 
<//>  is  the  maximum  stratospheric  heating  rate,  ti  gives  the  maximum  polar  night  jet  strength,  and  tp'h  the  gcopoten- 
tial  height  associated  with  planetary  wave  numbers  one  and  two  at  SO  km  and  S0°N. 
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